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CN ! Abstract 

An explicit solution of non-critical time-homogeneous branching processes is de- 
scribed. 
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\Q . 1 Introduction 
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Branching processes are widely used in high energy physics [1]. For example, 
the well known Furry- Yule and negative binomial distributions occur in simple 
branching processes with allowed transition 1 — > 2. The use of processes with 
higher order transitions 1 n with n > 2 is rare due to the absence of explicit 
solutions in terms of elementary functions. In this paper we describe the solu- 
tion for processes with higher order transitions using a recently found recursive 
procedure for the pure birth branching process[2]. The solution is based on 
the use of the Koenigs function[3] and the functional Schroder equation [4], 
sometimes called the Schroder-Koenigs equation (for a detailed description 
and extensive bibliography see [5,6]). In section 2 we describe the solution [2] 
for the pure birth branching process. In sections 3 and 4, respectively, the 
procedures for non-critical branching processes and branching processes with 
immigration are outlined. The results are summarized in the last section. 
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2 Solution for the general pure birth branching process 



A branching process with continuous evolution parameter t is determined by 
the rates a n for the transition ( "splitting" ) of one particle into n particles with 
all particles subsequently evolving independently. For a pure birth branching 
process «o = 0. The probability distribution p n (t) for the process having one 
particle at t = with p n (0) = 5\ n is a solution of the forward Kolmogorov 
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equation [7,8] 

dm dm 

-at = f(x) ^ 



for the probability generating function 

oo 

m(x,t) = J2Pn(t)x n (2) 

?1=0 



with 

oo 

f( x ) = a ^ n - ax ( 3 ) 

n=2 



and with a = J2 a n- The Taylor expansion of equation (1) leads to the following 
system of equations for the probabilities p n 



dpi m 

~dt = ~ aPl ' (4) 

= a 2 pi - 2ap 2 , (5) 

^ = a 3 pi + 2a 2 p 2 - 3ap 3 (6) 
at 

and for arbitrary n 

dpn \ ^ . /_•> 

^ = Z^^ a ™-i+iPi ~ • ( 7 ) 



A simple interpretation of the equation (7) is as follows: let us consider the 
state with n particles at the moment t. The change in this state is due to 
the arrival from states with multiplicity lower than n and to the departure to 
states with higher multiplicity. The arrival rate from the state with j particles 
is proportional to the number of particles j, to the transition rate (for one 
particle) to produce n — j new particles, i.e. a n _ J+ i, and to the population 
density in the state j, the sum in the equation (7) goes over all states below n. 
The departure rate is proportional to the total transition rate (for one particle) 
a, to the number of particles n in this state and to the population density, this 
explains second term in the equation (7). Formally this system of equations is 
valid for any initial condition. 
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Let us recall the recursive solution of the equations (4)-(7) given in [2]: the 
probability Pi(t) is pi = exp(— at) and 



Pn = J2 ^Pi (8) 
i=i 



with the following recursion for the coefficients iTj n 

n-j 

(n - j)ir jn = J2(n- l)k vr i(n _ . (9) 
i=i 

Here b\ = a\ + \/a is the relative probability to produce / new particles. The 
recursion starts from ir n = 1 and the coefficient ir nn can be found from the 
initial condition 

n-l 

TTnn = - 71 3™ ' ( 10 ) 

i=i 



For the case with N initial particles with p^(0) = Snu, the first N-l equations 
are automatically valid since 

^=^ = ...=P^! = (11) 



and the solution in this case has the following form: pffl = exp(-Nat) = p± 
and 

n 

pL n) = E *j?Vi (12) 

j=N 



with the same recursion as (9) for the coefficients 



(n - J)nff = n f(n - l)k . (13) 

1=1 



This recursion starts from = 1 and the coefficient irtfp can be found from 
the relation 

n—l 

(A ° - - E ■ (14) 
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j=N 
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One can calculate the coefficients ir^ using the concept of the Koenigs func- 
tion [2,5,6,9]. For the branching process starting from one particle at t — 
this function is defined as the limit 

tw \ i. m(x,nt) 

K(x) = hm — ^ — - . (15) 



K(x) has the following Taylor expansion: 



oo 



K(x) = KjX° = ]T ^3 X ° ■ ( 16 ) 

3=1 3=1 



For the branching process starting from N particles the Koenigs function is 
defined analogously: 

U'l J j=JV j=N 



The recursion (13) leads to the following recurrence for the coefficients N = 
1,2,...; j = N + l,N + 2,... : 



Let us denote n^+ n = t n (x) , then t (x) = 1 and t n (x) is given by the following 
recursion 

n 

nt n (x) = ^(x + n - l)b t t n -i(x) . (19) 
i=i 

It is evident from the equation (19) that the t n (x) is a polynomial of order n 
in x. The nf^ in terms of t„(x) is equal to tj_jv(iV). 

The remarkable property of the Koenigs function is that it satisfies the func- 
tional Schroder equation: 

K{m) = piK(x) . (20) 



It is convenient to introduce the function Q(x) the inverse of the Koenigs 
function. Then the equation (20) gives the functional relations m(x,t) = 
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Q(piK(x)) and m N (x,t) = Q N (p\K(x)) . The coefficients of the Taylor ex- 
pansion for Q N (x) = Qf^x^ can be found using the following relation (see 
Appendix B in [10] and references therein): 

Qf ] = ■ (21) 



In terms of t n (x) it gives 



Q? ] = jtj-NH) ■ (22) 



Finally, comparison of equation (12) with the Taylor expansion in x of the 
Q N (pi(K(x)) leads to the following expression for the coefficients nf^: 

vrg = Qf>«J?) = y/ ; vi j)t n ,!./) • (23) 



3 Solution for non-critical branching processes 



In this case the absorption coefficient «o is not equal to zero and the function 
f(x) (3) has an additional term ao(l-x). Let us denote by (5 the smallest pos- 
itive, different from unity root of the equation f(x) = 0. Branching processes 
with < f3 < 1 are called supercritical branching processes and processes with 
1 < (5 are called subcritical ones. 

The solution for the supercritical branching processes is obtained in the fol- 
lowing way. Let us transform m and x to m' and x' by the linear transform 

/ x-/3 , m-/3 . . 

X = —p ' m = — ' ^ 



this transform moves the root P to zero. Then the forward Kolmogorov equa- 
tion (1) for the process with absorption transforms to the form for the pure 
birth branching process. Therefore one can introduce the Koenigs function 
with 



K 
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where q\ = exp (— a't) . This leads to the expression for the m(x,t) 

m(x,t)=P + (l-P)Q(q 1 K(^y) (26) 

and to the infinite series in qi for the probabilities p n 

oo / ff\l-n 7| I 

^fc + li-fflEpL^Efttf . (27) 

The solution for the subcritical branching processes {(3 > 1) can be obtained 
in a similar way. In this case the linear transform should move 1 to zero and 
P to 1, i.e. x' — (x — I)/ {f3 — 1). This leads to the similar expressions 

m(x,t) = 1 + ((3 - l)Q^q lK (?0^j (28) 

and 

oo / 1 y— n l\ I 

The procedures described above are not suitable for critical branching pro- 
cesses with (3 — 1. In simplest case without higher order transitions the prob- 
ability generating function m(x, t) is a linear fraction with coefficients having 
linear dependence on t (instead of exp (—at) dependence for non-critical pro- 
cesses). In general case with higher order transitions the interplay between t 
and exp(— at) dependences leads to more complex equations. 



4 Solution for branching processes with immigration 



For the branching processes with immigration there is an additional external 
source of particles appearing in clusters of j particles with the differential rates 
(3j ( 12 13 j — b ). The generating function for the process starting with zero 
particles at t — can be written [11,12] as 

t 

M(x,t) = exp^y g(m(x,T))dT^j (30) 
o 
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with 

oo 

<7(s) = £/W-& , (31) 



i=i 



where m(x, r) is the solution for the underlying branching process without 
immigration. For the underlying pure birth branching process, equation (30) 
leads to the following expression 

M(x, t) = exp (-bt) exp (jr v n (t)x n ) (32) 

\i=i ' 



with 



•'n(t)=±±«®l-^ - (33) 

felj=1 J a 



Let us denote 



exp(J2v n x n ) = l + J^V n x n , (34) 

n=l ra=l 

then the final probability P n (t) = exp (— bt)V n (t). The coefficients V n (t) can 
be calculated using the recursive relation: 

n 

nV n = ^2jVjV n -j . (35) 

.7 = 1 



This relation is known in combinatorics and is used, for example, in the study 
of combinants[13-16]. 



5 Summary 



In this paper we have derived explicit expressions for the probability distri- 
butions in various branching processes. Although we have not derived closed 
expressions for the polynomials t n (x), the given recursions can serve as a cal- 
culational tool both in theoretical and experimental studies. 
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